Minimum-error discrimination between subsets of linearly dependent quantum states 
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A measurement strategy is developed for a new kind of hypothesis testing. It assigns, with mini- 
mum probability of error, the state of a quantum system to one or the other of two complementary 
subsets of a set of N given non-orthogonal quantum states occurring with given a priori probabili- 
ties. A general analytical solution is obtained for N states that are restricted to a two-dimensional 
subspace of the Hilbert space of the system. The result for the special case of three arbitrary 
but linearly dependent states is applied to a variety of sets of three states that are symmetric and 
equally probable. It is found that, in this case, the minimum error probability for distinguishing 
one of the states from the other two is only about half as large as the minimum error probability 
for distinguishing all three states individually. 

PACS numbers: 03.67-a, 03. 65. Ta, 42.50.-p 



I. INTRODUCTION AND BASIC EQUATIONS 

Due to thcir nonvanishing mutual overlaps, non- 
orthogonal quantum states cannot be perfectly distin- 
guished. Howevcr, stimulatcd by the ràpid developmcnts 
in quantum information theory 0], the question as to 
how to discriminate between non-othogonal states in an 
optimum way has gained renewed interest Q . In particu- 
lar, in quantum communication protocols several secure 
schemes have been suggested based on communicating 
via non-orthogonal quantum states. As a result, opti- 
mum discrimination between them became an inherent 
part of these schemes. For studying state discrimina- 
tion, it is assumed that a quantum system is prepared in 
one of the N pure states, \ipk), that belongs to a given 
set of non-orthogonal states, {\ipi}, IV^), ■ • ■ \iPn)} : an d 
that the a priori probabilities r]k for the preparation of 
either one of the states \ipk) are also known. In order 
to devise an optimum state-discriminating measurement, 
strategies have been developed with respect to various 
criteria |j, ||] . The earliest and simplest of these criteria 
is the requirement that the probability of getting a wrong 
result be as small as possible, with inconclusive results 
being forbidden and all states being individually distin- 
guishcd. A minimum-error strategy of this kind has been 
developed for the case when only two states are given 
and for specific N state problems || ^, [?], , includ- 
ing N symmetric and multiply symmetric [3| states. 
Recently the optimum strategy has also been found for 
three states exhibiting a mirror-symmetry || but still 
no exact solution has been known for N > 2 arbitrary 
states. Using the polarization states of a single photon, 
minimum-error discrimination has been experimentally 
realized for up to four symmetric non-orthogonal states 

In this paper we are concerned with a minimum-error 
strategy that involves N > 2 arbitrary linearly depen- 
dent quantum states, by considering the following prob- 



lem: We want to devise a measurement that allows us 
to decide, with the smallest possible error and without 
inconclusive answers, whether the actual state of the 
system belongs to the subset of states {IV'i)) ■ ■ ■ I^m)}, 
or to the complementary subset of the remaining states 
{|^>m+i)j ■ ■ ■ IV'Jv)} with M < N. For three given states 
the task reduces to distinguishing the state \ipi) from the 
set of states {(V^), IV^)} an d can be referred to as quan- 
tum state filtering with respect to the state This 
task has recently been investigated for the particular op- 
timization strategy that yields unambiguous discrimina- 
tion at the expense of allowing inconclusive results to 
oceur, the probability of which is minimized [ pd| . 

To treat our general minimum-error problem, we fol- 
low the Standard lines and introduce two positive Hermi- 
tian quantum detection operators, ITo and lli j^, ||]. We 
define the operator lli by the property that (ipk\H-i\Í>k) 
accounts for the probability to infer, from performing 
the measurement, the system to be in one of the states 
{\tpi), ■ ■ ■ iV'Af)}? if it has been prepared in the state \t[>k)- 
Obviously, this inference is incorrect if k > M. Simi- 
larly, given again the preparation of the state \ipk)i the 
quantity (-0^ |IIo I ^fc) denotes the probability for infer- 
ring the state of the system to belong to the subset of 
states {\ipM+i)i ■ ■ ■ IV'w)}; which is an erroneous result if 
k < M. Clearly, the relation 
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n, = i 



(i) 



has to be obeyed, where ï is the unit operator. From the 
definition of the detection operators it follows that the 
probability to get a correct result reads 
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p m{n) = Y,vM^k)+ w^fcinolV*). (2) 
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k=M+l 



In order to devise the desired minimum-error measure- 
ment scheme, we have to determine the particular detec- 
tion operators TIq and lli that maximize the right-hand 



2 



side of Eq. (||) under the constraint ([!]). In general, the 
error-minimizing optimization problem is a highly non- 
trivial task. 



II. SOLUTION IN TWO DIMENSIONS 

To enable simple analytical Solutions, we restrict our- 
selves to the case when the ./V linearly dependent states 
span only a two-dimensional Hilbert space. We note that 
for three linearly dependent states this is always the case. 
First, we show that in the two-dimensional case it is pos- 
sible to represent the two detection operators üi and n 
by two projection operators onto orthonormal states \fi) 
and respectively. To see this, we start from the ex- 
pression TIi = Ai|ui)(ui| + A 2 |v2)(^2|, with \v\) and \v 2 ) 
being the orthonormal eigenstates that belong to some 
non-negative eigenvalues Ai and A2. Expanding a par- 
ticular state \ip) as \ip) = cos/3|wi) -I- sin/3|^ 2 ), where a 
possible relative phase factor has been included into the 
definition of I «2), we arrive at 



(3) 



provided that we define \fi) = \/Aïl w i) i iyp^v^)- With 
the help of the relation 1 = + ^2) {^2 1, we obtain 

in the same way the representation 



(4) 



provided that \v) = \/l~~XÏ\vi} ± — Aal^)- Now 
we require that \(ii\tp k )\ 2 + KHV'fc)! 2 = 1 for an arbitrary 
state, \ip k ) — cosf3k\vi) + e' 7fc sinf3 k \v 2 ), which implies 
that |/í)(/í| + |^) {v\ = 1 has to be fulfilled. This only holds 
true when in the representations of \fi) and \v) opposite 
signs are chosen and when in addition A2 = 1 — Ai, lead- 
ing to the orthonormality conditions (/í|/í) = = 1 

and (fJ.\i > ) = 0. Therefore, in a two-dimensional Hilbert 
space the optimization problem posed by Eqs. (|ï|) and 
(|J) can be reduced to the problem of finding the specific 
normalized state that maximizes the expression 
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m (l-\(^ k )\ 2 ), (5) 



which follows when TEi = and IIo = ï — are 

substituted into Eq. (g). Comparing this to the spectral 
representation of the detection operators, introduced be- 
fore Eq. (||), we are led to identify |^í) with \vi) and \v) 
with \v 2 ) since the representation is unique. Then Ai = 1 
and A 2 = follows. Once the optimum detection state is 
known, the maximum achievable probability of correctly 
assigning a quantum state to one of the two subsets, as 
well as the two detection operators necessary to perform 
the optimized measurement, are uniquely determined. 

To solve the optimization problem, it is convenient to 
write the overlaps between the given states as 
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\A kl \e 1 ' 



(6) 



and to introduce the auxiliary state vector 
, , 1 



v/1-1^121 2 



(|V>2)-A 12 |^)). 



(7) 



For Eq. (||) to be vàlid, we have to assume that all N 
given states lie in a two-dimensional subspace, spanned 
by the states \ipi) and \ip 2 ), or j^i) and \v), respectively. 
Since (v\v) — 1 and (ipi\v) — 0, the states and \v) 
provide a suitable orthonormal basis for representing any 
state, iV'fc); as 



\ipk) = A lk \ipi) 



with 
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y/l-\A 12 \*y/ï=\A 



(8) 
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The last equation can be verified by calculating the over- 
lap (ip 2 \4'k}, taking into account that 72 = because 
of the specific definition of the state |t>). Similar ly, we 
represent the detection state, |/lí), as 



\fi) = cos i/? | Vi) + elx sinyj \v), 



(10) 



and obtain 



( M |V fc ) = Autcos^ + e^-^Vl- l^ifcpsin^. (11) 

Eq. ( |ï(ï| ) accounts for all possible states in the two- 
dimensional Hilbert space of interest provided that both 
tp and x are variables in the interval [0, 7r). The error- 
minimization problem is then reduced to finding those 
vàlues of ip and \ m Eq. (O) that maximize the proba- 
bility P M(N 1((p, X ) in Eq. (§). 

The solution to this optimization problem is straight- 
forward. We begin by inserting Eq. ( |ÏI| ) into Eq. (||) 
and, by making use of the fact that the a priori proba- 
bilities of the states fulfill the relation X^fcLi Vk = 1, we 
readily arrive at 

P M(N) = Í+i?cos(2^) + |Q|sim>)cos( X -XQ), (12) 
where R and Q are defined as 



R = 



and 
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(13) 
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(14) 



The conditions for an extremum, dP MÍ " N ^ / dp — and 
QP M ( N ) I = o, hold for p — tp c and \ — \ c , with 
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sin(2^ c ) = \Q\/y/W+W, cos(2^ e ) = R/y/& + \Q\*, 
and Xc = XQ: respectively. Note that cos(2</j e ) and R 
have the same sign while sin(2</? c ) is always positive. 
This choice of \ c and ip c corresponds to the maximum of 
P M W and, from Eq. @, we obtain 



p m(n) (^x c ) 



p M(N) _ i 
max 2 



VR 2 + \Q\ 2 - (15) 



The corresponding detection state, onto which a pro- 
jection has to be performed in a measurement scheme 
achieving the maximum probabüity, is determined by 
|/Xe) = cos(,5 |V>i) + e ÏXo sin<^ c |«). 

As applications of this general expression, we discuss 
two special cases. First, the solution can be cast to a 
considerably simpler form when the states are real. Real 
states have been considered before [[L2| in a different con- 
text. In this case the parameters of the optimum detec- 
tion state, can be calculated very easily. Both R and 
Q are real, yielding \c — if Q > or n if Q < 0. The 
maximum probabüity of determining correctly to which 
of the two complementary subsets a state belongs is given 
in this case by Eq. (|Ï5|) with 
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and 
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N 



fe=M+l 



(16) 
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where, in the last step, we made use of the relation re- 
sulting from Eq. ([}]) with jk = and all the overlaps are 
assumed real. 

As our second example, we consider the case of three 
arbitrary but linearly dependent states, N = 3. Choosing 
AI = 1 and taking 171 + 772 + 7/3 = 1 into account in Eq. 
(|Ï3|), we readily obtain 



R = --m\A 12 \ 2 -m\M3\ 2 . 



(18) 



The evaluation of \Q\ is greatly facihtated if we notice 
that the first sum on the r.h.s. of Eq. ( |Ï4| ) has only one 
term and this term vanishes. A straightforward evalu- 
ation of the remaining two terms from the second sum 
yields 

|Q| 2 = ^|A 12 | a (l-|^ia| 2 ) + ^|Aia| a (l-|Ai3| a ) 

+2ï72%(ReAi2A 23 A3i - \A 12 \ 2 \A 1Z \ 2 ). (19) 

For this case the parameters of the optimum detection 
state |/t) can be seen to be \c = XQ an d tan(2(^ c ) = 
\Q\/R, with \Q\ and R substituted from the above equa- 
tions. We do not give here a more explicit expression 
for Xc because it is slightly involved and enters only the 
detection states but not the final result for the maximum 
probabüity Inserting the above vàlues of |Q| and R into 



the general expression for the optimum probabüity finally 
gives 



pl(3) = 1 + 1 

max 2 2 



l-4E%(l-^)|(^i|V'fe)| 2 



k=2 



+8 r/2773 Re ((^1 ^2) (H^a) W'alV'i))] 



(20) 



This expression describes the maximum attainable prob- 
abüity of correctly distinguishing the state from the 
set of states {\tp2), \ipa)}- The minimum error probabüity 
then follows as P^rror = 1 — Pmax • As expected, the result 
is independent of the individual phase factors of the given 
states, and for 7/3 = it reduces to the pioneering for- 
mula jj] for minimum-error discrimination between only 
two non-orthogonal states. 



III. DISCUSSION 

With respect to possible applications, the question 
arises how the maximum probabüity for getting a correct 
result in quantum state fütering compares to the maxi- 
mum probabüity for correctly discriminating, by means 
of a different measurement strategy, between all the given 
states individually. In the following we shall explore this 
question for a variety of symmetric states. 

Let us investigate the set of three symmetric states 



\ipk) = cos/3|ui) +e' 3 



(fc-i) 



sin/3 |u 2 ), 



(21) 



with k = 1,2,3 and < (3 < 7r/4, which are assumed to 
oceur with equal a priori probabüity. Here \ui) and |it 2 ) 
denote any two orthonormal basis states. Obviously the 
states are linearly dependent and non-orthogonal. Duc 
to their symmetry, the mutual overlaps are equal and we 
get 4|A / t;| 2 = 4 — 3sin 2 (2/3) if k ^ l, where we again used 
the abbreviation Aki = (ipu^i)- Moreover, we obtain 
that 8Re(Ai2A 23 A3i) = 8 - 9sin 2 (2/3). By substituting 
these expressions into Eq. ( p0| ) and taking into account 
that rjk = 1/3, we find the minimum error probabüity for 
quantum state fütering with respect to the state \ipi), 



P 1{3) (8) - - 



v /l + 3sin 2 (2/3) 



(22) 



Because of the symmetry, the same expression holds for 
distinguishing any other state from the remaining two 
states. For comparison, we now consider individual dis- 
crimination between all three states. The general formula 
for minimum-error discrimination between N symmetric 
states, derived in Ref. H, has been recently applied by 
one of us |Ï3| to states of the form (^ï|), yielding the max- 
imum probabüity Pmax' 3 "* = |(| sin/3| + | cos/3|) 2 for cor- 
rectly distinguishing each state individually. From this 
result we obtain the minimum error probabüity 



P, 



(1,2,3) 



Error 



(/3) = l-P^· 3 ) = -[2-sin(2 / 3)]. 



(23) 
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The ratio Pj^^iP) / P^ r fo^ (P) ^ s f° un d to vary between 
0.5 for [3 = or n/4, and the maximum value 0.56 for 
j3 m w/12. When fi approaches zero, the physical diffcr- 
ence between the states vanishes and the respective mini- 
mum error probabilities, corresponding to random guess- 
ing, are twice as large as those for (3 = 7r/4, when both 
kinds of minimum error probabilities take their smallest 
possible vàlues. These vàlues are equal to 1/3 when all 
three states are discriminated individually, and to 1/6 
when only one of the states is distinguished. 

The same vàlues of the respective minimum error prob- 
abilities also result for the set of equally probable real 
symmetric states = \ui), [fa) = ~^(\ u i) + V^l^)), 
and ^3) = —l(\ui) — V3\u2)) which are known as the 
trine states jHj. For the case that \ui) and {112} refer 
to a single photon and represent horizontal and vertical 
linear polarization, respectively, these states have been 
used to verify experimentally the theoretical result 1/3 
for the minimum error probability in individual statc dis- 
crimination |]ïo) . On the other hand, from Eq. ( po| ) with 
r/k = 1/3 we easily find that the minimum error proba- 
bility for distinguishing the state alone is only 1/6. 
By using tan(2v3 ) = \Q\/R and Eq. |Ü]) with x = 0, 
the proper projection state, \p c ), is found to be \u\). 
Hence the corresponding quantum-state-filtering experi- 
ment for single photons could be performed with the hclp 
of a polarizing beam splitter that transmits the horizon- 
tal component and reflects the vertical one, or vice versa, 
as it is immediately expected in view of the symmetry of 
the problem. 

In conclusion, we remark that it is straightforward to 
generalize our bàsic equation (Q) in order to account for 
discrimination between more than two subsets. However, 
since the detection operators always have to resolve the 
idcntity, they cannot be represented by projection opera- 
tors onto orthogonal states if their number is larger than 
the dimensionality of the underlying Hilbert space. The 
measurement therefore would be a generalized |Ï4j mea- 
surement in this case. The same applies if the number of 



detection operators is smaller than the number of dimen- 
sions of the Hilbert space, as it happens if, e. g., Eq. (||) 
is applied to three linear ly independent states. Finally it 
is interesting to relate our results to the Helstrom bound 
Pe = — H w iPi™ w 2P2||)] for the minimum error prob- 
ability of discriminating between two density operators 
pi and P2 having the a priori probabilities wi and W2 , re- 
spectively. Here the symbol 1 1 • 1 1 denotes the trace norm 
||er|| = TrVcrto-. After inserting w 1 p 1 = Y,k=i Vk^k) {Í>k\ 
and W2P2 = J2k=M+i r lk\'4>k)('4>k\, the expressions ensu- 
ing from Pg for the cases we are interested in indeed 
confirm our results, without yielding the optimum detec- 
tion operators, however. 

To summarize, we derived the measurement strategy 
that minimizes the error probability for discriminating 
between two complementary subsets of a set of N non- 
orthogonal quantum states spanning a two-dimensional 
Hilbert space. The corresponding measurement is found 
to be a Standard von-Neumann measurement, projecting 
onto two orthonormal states that have been determined 
in the paper. Assuming arbitrary a priori probabilities of 
the N linearly dependent non-orthogonal states, we ob- 
tained a general analytical expression for the minimum 
error probability or, equivalently, for the maximum prob- 
ability of obtaining a correct result. As special cases of 
this general result, we gave explicit expressions for the 
case of N real states, Eq. ([ï^), and for three arbitrary 
states, Eq. (§(j). 



Acknowledgments 

J. B. wants to acknowledge the hospitality extended to 
him during his stay at the Humboldt-University in Bcrlin. 
The research of J. B. was also supported by the Office of 
Naval Research (Grant Number: N00014-92J-1233) and 
by a grant from PSC-CUNY. 



[1] 
[2] 

[3] 
[4] 
[5] 
[6] 

[7] 
[8] 



M. A. Nielsen and I.L. Chuang, Quantum Computation 
and Information (Cambridge University Press, 2000). 
For an overview on quantum state discrimination, see the 
review article A. C hefles, Contemp. Phys. 41, 401 (2000), 
quant-ph/00101l4 

S. M. Barnett, C. R. Gilson and M. Sasaki, J. Phys. A: 

Math. Gen. 34, 6755 (2001), |quant-ph/010702l| . 

C. W. Helstrom, Quantum Detection and Estimation 

Theory (Acadèmic Press, New York, 1976). 

H. P. Yuen. R. S. Kennedv. and M. Lax. IEEE Trans. 



[9] E.Andersson. S. M . Barnett, C. B. Gilson, and K. Hunter, 
|quant-ph/0201074 . 
[10] R. B. M. Clarke, V. M. Kendon, A. Chefles, S. M. Bar- 
nett, E. Riis, and M. Sasaki, Phys. Rev. A 64, 012303 
(2001). 

[11] Y. Sun, J. A. Bergou, and M. Hillery, |quant-ph/0112051 

(submitted to Phys. Rev. A). 
[12] M. Sasaki, S. M. Barnett, R. Josza, M. Osaki and O 
Hirota, Phys. Rev. A 59, 3325 (1999) 



[13] U. Herzog Fortschr. Phys. 49, 981 (2001), quant 



Inf. Theory bf IT-21, 125 (1975). 



M. Ban, K. Kurokawa, R. Momose and O. Hirota, Int. J. 
Theor. Phys. 55, 22 (1997). 

M. Sasaki, K. Kato, M. Izutsu, and O. Hirota, Phys. Rev. 
A 58, 146 (1998). 

S. M. Barnett, Phys. Rev. A 64, 030303 (2001). 



ph/0105139 



[14] K. Kraus, States, Effects and Operations: Fundamen- 
tal Notions of Quantum Theory (Springer-Verlag, Berlin, 
1983). 



